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EFFECT OF VORTICITY COHERENCE ON 
ENERGY-ENSTROPHY BOUNDS FOR THE 3D 
NAVIER-STOKES EQUATIONS 

R. DASCALIUC 1 , Z. GRUJIC 2 , AND M. S. JOLLY 3 


Abstract. Bounding curves in the energy,enstrophy-plane are derived for the 
3D Navier-Stokes equations under an assumption on coherence of the vortic- 
ity direction. The analysis in the critical case where the direction is Holder 
continuous with exponent r = 1/2 results in a curve with extraordinarily large 
maximal enstrophy (exponential in Grashof), in marked contrast to the sub- 
critical case, r > 1/2 (algebraic in Grashof). 


Dedicated, to Ciprian Foias in appreciation of his unbounded generosity 

1. Introduction 

Numerical simulations of turbulent flows (cf. [Q |2TJ [271, 3l|) reveal regions of 
intense vorticity dominated by coherent vortex structures; more specifically, vortex 
filaments. One of the imminent morphological signatures of the filamentary geom¬ 
etry is the local coherence of the vorticity direction ; it turns out that this property 
of turbulent flows leads to the geometric depletion of the nonlinearity. 

The pioneering work in this direction was made by Constantin. He obtained 
in [S] a singular integral representation of the stretching factor in the evolution 
of the vorticity magnitude featuring a geometric kernel that is depleted by local 
coherence of the vorticity direction, a purely geometric condition. This led to the 
first rigorous confirmation in J[9] of the local anisotropic dissipation in the 3D NSE: 
a theorem stating that as long as the vorticity direction is Lipschitz-coherent (in 
the regions of high vorticity), the L 2 -norm of the vorticity is controlled, and no 
finite time blow-up can occur. 

Subsequent work delved further into this geometric condition. The relaxation of 
the Lipschitz-coherence condition to a ^-Holder condition was made in [3], followed 
by a full spatiotemporal localization of the ^-Holder condition in [I]5|. A family 
of local, hybrid, geometric-analytic regularity criteria including a scaling invariant 
improvement of the ^-Holder condition was presented in ED]. Studies of coherence 
of the vorticity direction-type regularity criteria on bounded domains in the cases of 
no-stress and no-slip boundary conditions were presented in [3] and [2], respectively. 
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In this paper we derive a bounding curve for the weak attractor in the plane 
spanned by energy e and enstrophy £ under the ^-Holder condition in [3j. These 
curves suggest two scenarios in which the Taylor wave number kt can satisfy 

(1.1) KT:= (W) ** 

where k is the largest wavenumber in the force. One is where both energy and 
enstrophy are small; the other is where there are brief, intense bursts of large en¬ 
strophy. It is shown in mm that EU guarantees the energy cascade predicted 
by the Kolmogorov theory of turbulence [22] , The 1-Holder condition is essential 
for our analysis, of course. Yet, if the resulting curve turns out to be sharp, this con¬ 
dition still allows for excursions with extraordinarily large enstrophy; the maximum 
value achieved by this curve is £ ma , x = 0(exp(G 2 )), where G is the Grashof number. 
This would mean that if the ^-Holder condition were valid, solutions that are in 
fact regular might appear to blow-up in finite time under numerical simulation. 

The approach for the bounding curves follows that of Foias and Prodi (see pages 
59-61 in [14] as well as mm)- Each curve solves an ordinary differential equation 
in which e, £ are the independent and dependent variables respectively. The ODE is 
formed by the quotient of estimates for the growth rates de/dt and d£ / dt in a region 
of the e, £-plane where de/dt < 0. This technique was applied to the 2D NSE in both 
the energy,enstrophy-plane in m, as well as in the enstrophy,palinstrophy-plane in 
m, the latter being more relevant to the Kraichnan theory of 2D turbulence [23] . 
Curves which close were also found in this manner in m for the 3D Leray-a and 
the 3D Navier-Stokes-a models, whose solutions are global in time and approach 
those of the 3D NSE as the filter width parameter a —> 0 [7[ [32]. As expected, 
however, these bounding curves for the a-models blow up as a —> 0. 

For comparison, we include in this paper curves that somewhat limit solutions 
for the full NSE, i.e., without any conditions. Of course, using current techniques 
for estimating growth rates, those curves do not close. In [26j Doering uses the 
system of differential inequalities given by growth rate estimates to show that for 
initial data satisfying £q < z/ 4 /eo the enstrophy remains bounded. The maximal 
enstrophy growth rate for given enstrophy was studied numerically in [25]. 

In order to properly assess the significance of the ^-Holder condition in the 
realm of dynamics within the energy, enstrophy-plane, we also include the results 
describing the effects of the sub-critical coherence case where the Holder exponent 
satisfies l/2<r<l, as well as the effects of an example of the more traditional 
scaling-invariant conditions on the vorticity magnitude, namely, the small, uniform- 
in-time, L 2 -spatial integrability condition. It is worth observing that there is a 
qualitative jump between the maximal enstrophy allowed in the sub-critical and in 
the ^-Holder coherence scenarios; more precisely, from algebraic in G to exponential 
in G. Perhaps more intriguing is the comparison between the whole range of the 
Holder coherence conditions and the L 2 -condition. With respect to the scaling 
inherent to the 3D NSE, the L 2 -condition is scaling invariant, and the Holder 
coherence conditions are all sub-critical; hence-in this metric-the L 2 -condition is 
a weaker condition. In contrast, the L 2 -condition produces the most restrictive 
maximal curve in the energy, enstophy-plane, and is—in this sense-a dynamically 
stronger condition. 
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After some preliminaries in Section 0 we derive the bounding curves for the 
full NSE in Section [31 The precise formulation of the Holder condition is given 
in Section [31 The critical case, Holder exponent r = 1/2, is treated in Section [5] 
resulting in a bounding curve in terms of incomplete gamma functions, followed 
by sharp estimates for the maximum enstrophy value for this curve. Section [7] 
and Section [8] address the sub-critical Holder coherence and the scaling-invariant 
L? -spatial integrability cases, respectively. 


2. Preliminaries 


We consider the incompressible Navier-Stokes equations (NSE) in a bounded, 
connected, open set SlgR 3 

du 

(2.1) — — vAu + (u ■ V)it + Vp = F , x G 52 

(2.2) V • v = 0 , 


with either periodic or no-slip boundary conditions. The NSE can be written in 
functional form as 

du 

(2.3) ——b vAu + B(u,u) = f , uGH 

at 

where H is an appropriate (A 2 -based) Hilbert space, A = —PA is the Stokes 
operator, B(u,v) = P(u ■ V)v, f = PF, and P is the Helmholtz-Leray projector 
onto divergence-free functions. In this case, A is a positive definite self-adjoint 
operator with compact inverse, and we will denote by A(> 0) its smallest eigenvalue. 
To measure the size of the body force /, we introduce the dimensionless Grashof 
number 



where || • ||2 is the L 2 -norm on 52. A more detailed description of this functional 
setting for the NSE can be found in GDI- 

One of the main results of the (yet incomplete) NSE regularity theory is the 
existence of Leray-Hopf weak solutions to the initial value problem associated with 
& For an arbitrary T > 0 (fixed), any Leray-Hopf solution is 2?(A)-valued 
except, possibly, on a closed subset of [0,T] of 1/2-dimensional Hausdorff measure 
0. Thus it is regular on an open set which can then be written in the worst case as 
a union of countably many mutually disjoint open intervals, a set that is dense in 
[0, 00 ). See [TO], [21], and [29] for general background on the regularity theory of 
the NSE. 

We will use several well-known properties of the nonlinear term B 1 namely 

(2.4) B{u, v, v) = 0 , y u€ H , v GV = V(A 1/2 ) 
and 

(2.5) | (B(u, u), Au )| < c\\A l / 2 u\\l /2 \\Au\\l /2 , 

(see (6.18) and (9.22) in [TO])- Throughout the paper c will denote a generic, dimen¬ 
sionless constant, which will change from line to line, except when distinguished by 
a subscript. 
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Within each interval of regularity, we can take the L 2 -scalar product of (12.311 
with u, and then applying (|2.4I) . the Cauchy-Schwarz and Young inequalities, one 
has the standard differential relation for energy 

(2.6) + 2v\\A 1 / 2 u\\l < 2||/|| 2 ||u|| 2 = 2^ 2 A 3 / 4 Ge 1 / 2 < ^A^G 2 + v\\\u\\ 2 , 


which, together with the Poincare and Gronwall inequalities, implies that for regular 
solutions 


limsup ||u|| 2 < 

t—*oo 



v 2 G 2 
A 1 / 2 ' 


In fact this property holds for general Leray-Hopf solutions. 

More precisely, the long-time behavior of (12.311 is encoded in the weak attractor, 
Aw (first introduced in [18]), the subset of Leray-Hopf solutions that are globally 
bounded in H - forward and backward in time. This set weakly attracts all Leray- 
Hopf solutions as t —> oo. By using the Leray-Hopf energy inequality one can 
establish that A w is contained in a ball in H centered at 0 with radius vG / A 1 / 4 , 
i.e., 


(2.7) 


Aw C H(0, i/G/A 1 / 4 ) . 


In particular B{ 0, uG/A 1 / 4 ), is a (forward in time) invariant set for the Leray-Hopf 
solutions. Therefore, in what follows our study of energy-enstrophy relations will 
be focused exclusively on this ball. 

We conclude by noting that the above theory of Leray-Hopf solutions will be less 
relevant in the no-slip case studied in Sections m since the vorticity coherence 
assumption (|4.4[) combined with EU> imply that the solution is in fact regular, 
and therefore estimates like (12.611 will hold for all t. 


3. Bounds for the full 3D Navier-Stokes equations 

While an overall bound on enstrophy for the 3D NSE remains an open question, 
the partial regularity results described above will allow us to derive certain curves 
in the energy,enstrophy-plane which limit the long time behavior of Leray-Hopf 
solutions, without any additional regularity assumptions. 

In this section we consider periodic boundary conditions on H = [0, L] 3 , take the 
phase space H for (12.31) to be the closure in L 2 (fi) 3 of all divergence-free, mean zero 
trigonometric polynomials, and denote the energy, enstrophy and palinstrophy by 

e=IMli, E = ||A 1/2 u||2 , P = ||-4w||l ■ 

As noted before, we will concentrate on the globally bounded solutions, for which 
(12.71) holds, i.e: 


(3.1) 


e < e 0 



v 2 G 2 

A 1 / 2 


Now on each interval of regularity, take the L 2 -scalar product of (HOD with Au, 
apply (12.51) , and then proceed as in (12.61) to find 


1 dE 

2 t it 


v P 


||/|| 2 P 1/2 + cE 3/4 P 3/4 



V 


vP cE 3 j/P 

— - 1 - 5 - 1 - T~ 

4 j/ 3 4 


< 
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Again by the Cauchy-Schwarz inequality we have E 2 < eP, so that 

dE < 2v 3 \ 3 ' 2 G 2 + ^ - uP 


dt 


p3 p2 

< 2u 3 A 3 ' 2 G 2 + - v— 


(3.2) 

Note that dE/dt < 0 provided 

(3 ' 3) 6 - ^ (E) = 2u 6 A 3 / 2 G 2 + Cl E 3 

and that the nonzero critical number for $ is 


z 4 E 2 


(3.4) 


Ei = 


Cl 


1/3 


For some 77 > 1 consider the region where we have both 
(3.5) 
and 


v 2 \ 1 / 2 G 2 / 3 . 

:re we hr 
E > rji/X 3/A Ge 1/2 

> 2v 3 \ 3 / 2 G 2 , 


ciE 3 


or equivalently 
(3.6) 


E > E = 


\ i/3 

- j zz 2 A 1 / 2 G 2 / 3 = 2 - 1 / 3 E 1 . 

: have 

< 2cii/“ 3 E 3 - 77 u 2 A 3 / 4 Ge- 1/2 E a.e. t € [0, 00 ) 
ves 

< — 2(77 - l)i/ 2 A 3 / 4 Ge 1/2 a.e. t £ [0 , OO) . 


Cl, 


(3.8) 


Using (13.51) and (13.61) in (13. 21) . we have 
„ dE 

< 3J > Tt 

while (13.51) in (EH gives 
de 
dt 

To find a bounding curve we proceed as in HUE]. It is constructed as a 
solution to an ODE in the variables e, E. in decreasing e. As long as (13.51) holds so 
that the upper bound (13.81) is not positive, and as long as the bound in right hand 
side of (13.71) is nonnegative, the slope of this curve is determined by quotient of the 
two bounds. The bounding curve is then the solution to the Bernoulli ODE 

dE _ 2cn/- 3 E 3 -7 ? u 2 A 3 / 4 Ge- 1 / 2 E 
de 


—2(77 — l)i/ 2 A 3 / 4 Ge 4 / 2 


(3.9) 


1 77 E 


Cl 


(77 — l)u 5 A 3 / 4 G 


=1/2 


2 77 — 1 e 
whose solution can be written as 

(3.10) E = *(e) = {(jY Eo 2 + P (e 1 / 2 - e"“eJ /2+ “) } 

where 


-1/2 




77 -I 


and p = 


4ci 


(377 - 1 )tAA 3 / 4 g ’ 
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The function $ has a vertical asymptote at e = e*, where 


a+l/2 _ a+1/2 

*=* — e n 


*e§- 


We assume that (13.61) holds for initial data (e 0 , E 0 ) on the parabola 
(3.11) E = ^A 3 / 4 Ge 1/2 , 

so that in terms of 77 , this asymptote is at 


e*(??) = e 0 1 - 


2 (^- 1 ) 
377 — 1 \ 3 "- 1 


4ci?7 2 G 4 / 

We next solve (13.91) with initial data (ei, Ei) given by the maximum of dr in (13.31) . 
where Ei is as in (13.41) and 

2/3 

ei = 


I ^ „ 2 A -v 2 g - 2 / 3 . 


The solution curve with this initial data would intersect the parabola (13.111) at the 
value of e satisfying 

(3.12) e 1/2+ “ — 7 e “ -1 + 5 = 0 

where 

1 , c e i l/2+a 

&nd S = JE 2 _ 61 


7 = 


/3r) 2 i/ 2 \ 3 / 2 G 2 

Substituting the values for ei, Ei, one can show that 5 > 0 provided 

, 4 Cl /4\ 5/6 
V< + 3 V6\cJ ■ 

Dropping the 7 -term in (13.121) . we obtain as a lower bound for its root 

e 2 > tfdWA-^G 2 / 3 ■ 


For large G, we then have that > §, where (e, E) is the point of intersection of 
the line E = JE and the parabola in (13.111) . A qualitative sketch of the shaded region 
in which the weak attractor must lie is given in Figure [T| It is not surprising that 
this region is not closed, since global regularity is still an open question. As far 
as these estimates show, the sufficient condition (large kt, see ( 11 . 11 ) ) for an energy 
cascade can be met either by the solution spending most of its time near zero, or 
with repeated excursions into the intermittent region. 


4. Vortex stretching depletion 

We now show how a geometric condition of the direction of vorticity allows us to 
close the region for the weak attractor. For convenience in the analysis to follow, 
we now consider no-slip boundary conditions 

(4.1) u{x) = 0 , ieF = dfl . 

Taking the curl of (12.11) gives the equation for the vorticity u> = curl u 
dix 

(4.2) —-pAw + [u ■ V)w = (to ■ V)u + curl F , x £ Q 
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E 


Towards singularity 


= 0(G 2 / 3 ) 


0(C 



ei = 0(G _2/3 ) 


0(G 2 ) 


e = 4>(E) 


e 


Figure 1. Bounds for the full 3D NSE. Region I is recurrent. If 
the solution blows up, it must enter region II. If it enters regions 
III, IV it must enter region I. In region IV the enstrophy is always 
decreasing, while in regions II,III, IV, the energy is always decreas¬ 
ing. 

The evolution of the enstrophy is found by taking the scalar product of (14.21) with 
w in L 2 (Q) 



(4.3) 


The obstacle to proving global existence of solutions is the vortex stretching term 
((w • V)u,cu). One way to mitigate its effect is to impose a geometric condition on 
it in physical space. 

This amounts to a coherence assumption on the direction of the vorticity of the 
form 


| sin0(w(x), uj(y))\ < c\ r/2 \x - y\ r 


(4.4) 


where 0 is the angle between the two vectors, and | • | denotes either || • || R 3, or 
absolute value. 

5. The critical case 

Following the work of Constantin and Fefferman in jft[ S3, H. Beirao da Veiga 
shows in [5] that condition (14.41) with r = 1/2, along with the additional assumption 


(5.1) 
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for some constant <5 < 1 and some ip G L 1 (0,T), is sufficient for regularity of the 
solution to the NSE on the interval [0,T]. To do so he splits the vorticity as 


uj{x) = u/ 1 ) ( x) + w (2) (x) 


ui^(x) if |w(cc)| < 
uj^(x) if |w(a;)| > y\v 


for some y > 0. This leads to an expansion of the trilinear term into eight sum¬ 
mands: 


((u> ■ V)u ■ u>)(x) 


2 

y ] fca 0 -y(x) 

a,/3,7=1 


2 

E (c 

e jk i u} 

a,/3,7=1 


{x)u\^ (x) 


d 2 G(a 


dykdx 


V ), ,(7) 
OJ, 


(y) dy 


where 


e jkl 


! 1 (i,j, k) even permutation 

—1 odd permutation , and Q{x,y) 

0 i, j, k all equal 


1 

47 r|a; — y\ 


+ g(x,y) 


for some regular function g. The only summand that uses (14.41) is YZiti- Since 
we will modify slightly Beirao da Veiga’s estimate, we recall part of his argument. 
First he shows that 


|A^ 222 (a;)| < c 


M 2 >Qr)|V 2 >(;/)| 

\x - y| 3 


sin 6 (x 1 y)\ dy . 


Then, writing / for the Riesz potential 

I{x) = f \ui(y)\- - —, 

Jn \x — y \ 5 ' 2 

using (14.41) with r = 1/2, and Holder’s inequality, he finds that 


(5.2) f \K. 222 (x)\ dx < cA 1/4 f \u(x)\ 2 I(x) dx < cA 1/4 ||a7|| 2 ||u2|| 6 ||i'||3 

Jo, Jq 

As in [2], we apply the Hardy-Littlewood-Sobolev inequality [25] ||J|| 3 < c||w ||2 
and the Sobolev embedding Hq(SY) C L 6 (fl), but then use the Poincare inequality 
before applying Young’s inequality 


(5.3) 


/ \IC 222 {x)\ dx < cA 1/4 ||w|| 2 (A 1/2 ||a ;||2 + ||Vw|| 2 )||w||: 
In 

< cA 1/4 ||Va;|| 2 ||^||| 

< ei'||Vcj||f + cA 1 / 2 UiliU-i. . 


ev 


Otherwise, except for the inclusion of dimensional factors involving v and A, we 
treat the bounds for the remaining summands as in [2]. They amount to the three 
cases 


| f IC a 02 (x) dx\ < cyv\\\uj\\l , for (a,/3) ^ (2,2) 
Jn 
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and the remaining four 

| [ /C a /3i{x) dx\ < £is\\Vuj\\l + cs~ 3/5 (lj,\) 4/5 v 1/5 \\uj\\2 /5 \\u}\\1 . 

Jn 

For the remainder of the paper we take £ = ||w||| and V = ||\7co||| as enstrophy 
and palinstrophy respectively. Assuming that (Id.41) . (15.11) hold, we then have that 


1 dS- 

+ uV < ( 2 e + 5)vV + C 2 

( 5 . 4 ) 2 dt 


„l/5 X 1 / 2 ' 

+ V’) + i3/s(M) 4 ' /5 ^ 2 ' /5 H- £ 


ev 


11 curl F|| 2 f 1/2 . 


The point in applying the Poincare inequality in (15.31) is that otherwise the sum¬ 
mand A 3 / 4 !? 1 / 2 would appear in the brackets. To further simplify, we note that the 
summand involving £ 2 / 5 will dominate both the first summand in brackets as well 
as the term from the force provided 


(5.5) 

£ > £ m in = max 


£ 3 / 2 \ 1 / 2 lS 2 

M 2 


O' 


3 ) 5/2 , 


=-2/3 


C2 0/9 (mA) 8 / 9 U 2 / 9 


||curl U || 2 ° /9 


where we have tacitly assumed that ip £ L°°(0,T). Now choose e small enough to 
satisfy 


(5.6) 


2s H- S — p <C 1;. 


Note that 

£ = IMI 2 < 2||^4 1 / 2 tz||| < 2 ||u|| 2 ||Am || 2 = 2 || U || 2 ||A «|| 2 

= 2 1|"Lt|| 2 1|curl w || 2 < 2||u|| 2 V2||Vw || 2 = 2 3 / 2 e 1 / 2 P 1/2 


hence 


(5.7) -8P < -£ 2 /e . 

Gathering terms, we find that 

(5 8 ) —<( 02X1/2 - £ 2 , 6 c 2 (/xA) 4 / 5 u 1/5 f . 7 /5 

dt ~ \ ev 4e ) £ 3 / 5 

We now assume, in addition to (15.61) . that 
(5.9) v£ > qf\\ 2 e 1/2 . 


Using this in (12.61) . along with £ < 2||A 1 / 2 it||| we Hnd that 


(5.10) 


de 

— < 

dt 



As in Section [3] we assume (12.71) holds, and so, by combining (15.101) with (15.81) 
as indicated in Figure E3a) we can construct a bounding curve for A w by solving 


(5.11) 

(5.12) 

(5.13) 


d£ = (1-p _ 2c 2 A 4 / 2 \ _ 12c 2 (/xA ) 4/5 2/5 

de \ 2e ev 2 ) £3/5 j,4/5 

v 2 G 2 

£ (e 0 ) = £0, where e 0 := and 

£ 0 := max |4 ||/|| 2 e^ 2 , £ min } = max |4u 2 A 1 / 2 G 2 , £ min } 
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provided the right hand side in (15.811 is nonnegative, i.e, (e, £) is to the right of the 
curve 

I 5 / 3 ,,2 


(5.14) £ = e^y/^/V 


6 e 


e a - e 


e < e a = 


e(l — p) v 2 
4 c 2 A 1 / 2 


Notice that if G is big enough, e 0 = vG 2 j A 1 / 2 > e a , and so the solution £ = </?i(e) 
to (15.11115.1^1) is indeed to the right of the curve in (15.141) over an interval (e max , eo]- 
In fact, tpi is a decreasing curve on [e max ,eo] with global maximum at (e max , £ m ax) 
- the intersection with the curve in (15.141) . 

Recognizing (15.111) as a Bernoulli equation, we set 


f = £ 3/5 

to obtain 

(5 - 15 > I 

whose solution is 
(5.16) 


3(1 ~p) 
10 


b = 


6 c 2 AV 2 
5 £V 2 ’ 


and C = 


36c 2 (mA) 4 / 5 

5 e 3 / 5 jy 4 /5 


+ (“ i + b ) Z = ~ C , = (2^ 3 A 1/2 C 2 ) 3/5 . 


£( e ) = e“e 


a be 


eo'“e be °eo - C / s~ a e os ds 


—aJ)s 


To evaluate the integral we use an expansion of the incomplete gamma function 

(-*)" 


(5.17) 7 (a, 2 :)= / t“ : e T dr = z a ^ 

Jo 

so that with a = 1 — a 


n—0 


n!(a + n) 


/>e p—be 

/ s~ a e bs ds = (- 6) a - 4 / T~ a e~ T dr 
J e0 J—be 0 


= (-&) a [7(a,-&e) - 7 ( 0 , -fee 0 )] . 


The solution to (15.111) can thus be written as 

5/3 


To find a bounding curve for e < e max under the assumptions (15.51) and & 
observe that since the right hand side in (15.81) is negative, we need a lower bound 
on de/dt. From Theorem 1.1 in [ 6 ] (henceforth, £7 is convex) there exists a constant 
cq > 1 depending on the domain, such that E < cq,£. Using this in 

-2II/U2NI2 < ^ + MA 1/2 u\\i < J + 2 "cn£ , 

together with (15.91) . we may now write 

(5.18) Co = 1 + 4 cq , 


£ = tp i(e) = { e a 


„ —be 


e~ a e be ° 


tio-C 


\ ' 


(fee) 


^ n\{a 

71=0 ' 


- e n 


(be 0 y 


OO 

y 

n\(a - 

n =0 ^ 


n) 


and combine with (15.81) to arrive at an initial-value problem similar to (15.1UI5.12)) . 
but with the right-hand side of (15.111) divided by Co, and the initial condition 
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Consequently we obtain a bounding curve £ = 932 (e) = £ 5 / 3 (e), where £ is as in 
(15.161) , but with a , b, C, e 0 and £ 0 replaced by a/C n , b/C n , C/Cq, e max and £max- 
It is easy to see that 93 2 is an increasing, concave curve on (0,e max ), and that 
ip 2 (e) \ 0 as e \ 0. It follows that on (0, e max ] the curve £ = 93 2 (e) will stay to the 
left of that in (15.141) : there will exist e m i n G (0,e max ) for such that 93 2 (e m i n ) = £ m in 
and 9 ^ 2 ( 6 ) £min for all e G (e m i n , e max ]. 

In order to show that the assumption (15.91) is satisfied on [e m ; n , e max ], note that 
from (15.161) . with the constants modified by Cn and with arbitrary eo, 0 < e < eo, 


£(e) > 


„ a £( e o) 


Thus, if we start above the parabola (15. 91) . i.e. with £ 0 = e 0 /3 > 4||/|| 2 eo /2 /^, and 
to the left of the curve in (15.141) . then the corresponding curve ip 2 will satisfy 


v£ > v£ 0 



> 4||/|| 2 eg 



4||/|| 2 ei (^)* ^ >4||/|| 2 ei 


and consequently, the curve £ = 932 (e) will satisfy (15.91) for e G [0,e max ]. Therefore, 
£ = 9 ? 2 (e) will be a bounding curve for A w for e G [e m i n ,e max ]. 

On e G [0, e m i n ], the assumption (15.51) does not hold, and thus the construction of 
bounding curves depends on whether v\(fj, + ip)£ dominates 11curl F \\2 £ 1 ^ 2 in US2D- 
This potentially would introduce two more cases. However, in the case 11curl F || 2 
is big enough, e.g. 

\5/4 2 

(5.19) 11curl F || 2 > c 2 £ 3/4 - (p + |fylU) 9/4 , 

/ 1 

the condition £ < f max would yield ||curl F^f 1 / 2 as the dominating term, and so 
instead of (15.81) we will obtain 


(5.20) 


d£ < / c 2 A 1 / 2 
dt ~ \ £v 


1/(1 4e P) )g 2 + 6||curlF|| 2 g 1 / 2 , 


which together with (15.181) yields the following equations for the bounding curve 
on [0, e m i n ]: 


(5.21) 

(5.22) 


d£ 

de 


fl-P 


\2Cne Cqsv 


2 ^ 1/2 ) £ - 


12 


11curl F || 2 £ -1/2 


£(e m in) — 9 5 2(6min) — £n 


Solving this Bernoulli-type equation we obtain 

/ e T/e\“ 

x{e) = x 0 I ) exp (/3(e min - e)) + 7 / (jJ exp (/3(t - e)) dt, 

e 

„ 3c 2 X 1/2 18 


= £ s,a , X„ = ( £ (.„,„))^. « = |i^ 

Note that the above equations imply 


Cq.£V 2 ’ 


7 = 




||curl F || 2 


x(e) > xo 
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and therefore, since (15.91) is satisfied by the initial condition (15.221) . meaning 

Xo /3 = £(e min ) > 4 ||/|| 2 e^/i/, 

we have for eg [0, e m i n ] 


£ = (x(e)) 2/3 > Xq 73 


^mir 


2/3 

- Z 0 _„l /2 
1/2 C 


(e^h,)i-|° £ 4||/|hel/> _ 


i.e. the assumption (15.91) holds on e S [0, e m j n ], which also means that £ will stay 
to the left of the curve in (15.141) . Therefore, 


£ = ip 3 (e) = (x(e)) 2/3 


is a bounding curve for A w for e e (0,e m ; n ]. It is easy to see that ip 3 is also 
increasing, concave and approaches zero as e \ 0. In fact, 

(5.23) /53(e) = O ^e 2 ^^ as e \ 0 . 

We denote the complete bounding curve by 


^(e), 

® ^ [^maxj ®o] 

</> 2 (e), 

® ^ [^min?^max. 

^3(e), 

G G (0, 6 m i n ) . 


Since by Poincare inequality Aoe < E, and by [ 6 ] E < cq£ , we have the following 
bounding region for A w in (e, £)-plane: 

A e < £ < <p{e). 

with A = Xo/cq. 

Numerical solutions of the bounding curve ip are plotted in Figure 2(b). One is 
made by using a 4th order Runge-Kutta method to solve the ODE, the other two by 
truncating the series approximation for the incomplete gamma function. Though 
the Grashof number is quite small (G = 2), the maximum value £ ma x is at least 
G(10 35 ). Rigorous upper and lower bounds for £ max in the next section confirm 
the dramatic rise in this bounding curve. 

In particular, estimates of Section [G] show that as the Grashof number grows, 
£m a x ~ exp(c 3 G 2 ), e a ~ e max ~ G°, and therefore, we can deduce from (15.161) 
that e m i n ~ exp(—C 4 G 2 ), for appropriate G-independent constants C3, C 4 Thus as 
G — > 00 , the possible big values for enstrophy will be constrained in a narrow 
(compared to G) strip of energy values located in the proximity of e = 0. 


6. Estimates for the maximum value of £ = ip(e) 

For simplicity we examine the case v = A = 1, so that the bounding curve 
£ = </?i (e) in the critical case is the solution to 


d£ _ f 2c 2 1 - p \ _ 12c 2 p 2/5 
de \ e 2 e / e 3 / 5 ’ 

£(e 0 )=£ 0 = 4||/|| 2 e 3 / 2 =4G 2 . 













(5.8) 
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Figure 2. (a) Schematic plot indicating which bounds on time 
derivatives are used to determine the ODEs for <p i and </? 2 , (b) 
Numerically generated bounding curve for u = A = p = 1, e = .2, 
5 = .5, C 2 = G = 2, by Runge-Kutta (RK) and by truncating 
the series approximation for the incomplete gamma function to N 
terms. Estimates for £ max are in Section [6] 


6.1. A lower bound: 

Clearly, 


d£ 
de ~ 


2 c 2 
£ 


1 ~ P 
2 e 


£. 


Therefore, on [e max ,eo] the curve £ = ip(e) lies above the solution to 


( 6 . 2 ) 


d£ 

de 


2 c 2 1 — p 

e 2e 




e 0 = G 2 , £ (e 0 ) = £ 0 = 4G 2 , 


which is 


eo 


2c 2 


£(e) = £q ( — ) exp I —(e 0 — e) ) = 4G 1+P e 2 exp I —(G 2 — e 


1 —p ( 2c 2 


It is immediate from the ODE in (16.21) that the maximum for £ is achieved when 

e(l - p) 


£ — ^crit — 


4c 2 


and therefore the maximal value of the bounding curve £ = ip(e) satisfies: 


(6.3) 


V _____ ( 2c 2 


> 4G 1+P e cr ; t exp 


-(G 2 — e C rit) ) • 
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6.2. An upper bound: 


Looking back at the original equations (imp . we notice that if for p > 0 
( 6 - 4 ) ^S>^£ 2/5 , or, equivalently p > ’ 

then 

— > _ _ 1 ~ P \ £ _ 021 £ = _ ( ( 2 + ^) c 2 _ 1 ~ P \ £ 

de ~ \ e 2 e ) e \ e 2 e J ' 

Note if (16.41) holds for £ = £q, then it holds for the bounding curve £ = </(e) over 
e € [e max , eo]. Thus, for r\ big enough so that (16.41) holds at £ = So, on the interval 
[e mil x, eo] the bounding curve £ = <p(e) will lie below the solution to the following 
initial value problem: 

d£_ _ _ ( (2 + rj)c 2 _ 1 - p \ ^ 

(6.5) de V 6 2e / 

e 0 = G 2 , £ (e 0 ) = £o = 4G 2 . 


Obviously, this is exactly (16.21) with 2c2 replaced by (2 + 77 ) 02 - Therefore, 
(6.6) f max < 4 G 1 +p eJ^ ex P ^ + ^ C2 (G 2 - e crit )^ , 


where 


^crit — 


e(l - p) 


2(2 + rf)c2 

Note that e cia t ft 6 max ^ Ccrit, so e max ~ £. 

Taking £ = .2, fi = 1 G = C 2 = 2 in (|6.4[) we can choose rj = 4.33 to evaluate the 
expressions in m and (16.61) and find that 

5.83 • 10 35 < £ max < 9.13 • 10 110 , 


which is consistent with the plot in Figure [U(b). Naturally, since G was chosen so 
small, the upper bound on E max is much less sharp than the lower bound. We can 
make the upper bound sharper if we set £(eo) = 39311.12 O 4G 2 in (16.51) . Then 
we can choose p = 0.04 resulting in the upper bound E max < 6.99 • 10 39 , which is 
much closer to the lower bound. Note that as G increases, so does So, and thus the 
smaller we may take 77 to sharpen the upper bound. 


7. The sub-critical case 1/2 < r < 1 
Instead of (15.21) we have 


(7T) |/C 222 (a:)| dx < c \ r/2 |w(x )] 2 ^ |J^| 3 _ r d vj dx 

(7.2) < cA^HwHl/allfe* MIU , 

where 

h(x) = ' 

To estimate the convolution, we apply the Holder inequality for convolutions ac¬ 
cording to 


1 + — = —|— where p 

4 P q 


3 12 

3 — r ’ ^ 3 + 4r 
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followed by interpolation to obtain 

l-{~4r 3 — 4 r 

II h* |w||| 4 < c||/i||p )tu ||a;||q < c||w|| g < c||w||~||Vw||~ 


Using this and the interpolation estimate 

IMli/3 — c ll c °’ II2 /4 II Vco || 2 /4 

into (EH), we have 

P 3 _j _2 r 3 2 r 

I \IC 2 22(x)\ dx < cA r ' /2 ||w||^ _ ||Vu;|| 2 ^ 


< eu|| Vco||| + c 


A 2 


(ei/) 3_2r 

This gives a bound just as in (15.41) . except that 


6+4r 

l 2 1+2r 


A 1 / 2 


£ is replaced by 


A 2 


£T+2F 


EV ~ ~ \(£l>) 3 ~ 2r 

Thus, under a threshold condition 

(7.3) £>£, 

where £_ is defined similarly to £ m \ n in (15.51) so that the summand with the 2 
power dominates the others, we have 

(7.4) 
where 


l+2r 


d£ Z/(l — o) _ 9 C 3 + 2r 

— <--- —£ 2 + 12 —£ t +37 

dt 4e 2 


C= -6c 
v 


A 2 


[£v) 3 ~ 2r t 

We then combine (17.51) with (15.101) to find a bounding curve solving 

d£ 
de 


(7.5) 


1 P.g _ C£^ 


2 e 


over the interval [e, eo] where (17.41) is nonnegative. The solution to (17.51) is 


(7.6) 

where 


£ =«', = 1UJ £S+ ~ 

2 r- 1 


aC 


ctcr 


and a = 


_ CKCT — 1 


1 ~ P 


— e 


l/o 


1 + 2r 2 

From EH, we have that the maximum value, £ = ^i(e) is achieved where the 
bounding curve (1731) intersects the curve a£ a = Ce. An elementary calculation 
gives, for appropriate G -independent constants C\ and C 2 , 

£= [C 1 (G 2 ) a - aa + C 2 (G 2 ) 1 - aa ]^ J ^ = 0(G 2/a ) 


as G —> oo. 

For e < e, as long as (17.31) and (15.91) hold, we combine (17.41) with (15.181) to obtain 
a bounding curve £ = </> 2 (e) where <p 2 is the same as (j> in (17.61) except a, C, eo and 
£o are replaced by a/Cn, C/Cn, e and £. As in Section [H we can show that </> 2 
is an increasing concave curve that converges to zero at the origin and stays above 
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the parabola £ = 4 ||/||2 e 1 / 2 / over the interval [0, e] . Therefore 02 will intersect 
the line £ = £_ at some point e £ (0, e), and so 


£ = </> 2 (e) 


is a bounding curve for A w on [e, e]. 

The bounding curve for 0 < e < e is obtained analogously to <^3 in Section [5] 
In fact if £ < £_ and ||curl F ||2 is bigger than an appropriate Grashof-independent 
threshold similar to (15.1911 . then the term ||curl will dominate, resulting 

in the following equivalent of (15.21115.22)1 : 



(7.7) 

(7.8) 

This yields the curve 


£(e) = 02 (e). 



(7.9) 


with /3 = (1 — p)/(2Cfi). As in Section 0 03(e) is increasing on (0, e], convergent 
to zero at 0, and is above the parabola £ = 4 ||/|| 2 e 1 ! 2 /v. Also, 



(7.10) 


as e \ 0 , 


i.e. the same rate as in critical case - see (15.231) . 

Thus we obtain the following bounding region for A w in (e,£)-plane: 

! 0 i(e), ee[e,e 0 ] 

02 (e), e e [e, e) 

03 (e), ee( 0 ,e). 

These bounds are illustrated in Figure [3] 


e 


40 



35 


30 


25 


20 


15 


10 


e 


0.5 


1.5 


2.5 


3.5 


Figure 3. Bounding curves for subcritical case, r = .51, G = 2, 
C = 1. Region I is recurrent. In region III the enstrophy is always 
decreasing, while in regions II and III the energy is always decreas¬ 
ing. 
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8. Bounds from scaling-invariant regularity condition 

A natural question to pose at this point-for the sake of comparison-is what are 
the effects of the classical, analytical scaling-invariant regularity criteria on the 
vorticity magnitude on the dynamics considered in the energy-enstrophy plane. As 
an illustration, we present the case of the small-La spatial integrability condition, 

sup ||w(t)||i < e 0 , 

*€(0,T) 

for a sufficiently small constant co¬ 
in this case, the estimate on the vortex-stretching term reads as follows. 


\ui ■ Vw ■ u>\ dx < ||w|| 3 ||Vll|| 6 ||w ||6 


< C n M,Ml < on |M| # ( ||Vco||| + fcnIMI 


< eocn ^||Vw||2 + MMI2 

< coc'qu (II Vco||| + A||co|||) . 

Using this along with (15.11) . (15.71) yields 

% < - j(l - eocn - 5)— + 2uA(i/> + e 0 c' n )£ + 2||curl F|| 2 £ 1/2 . 
at 4 e 


Now assume 


so that 


£ > 


d£ 


||curl U|| 2 \ 2 

+ e 0 c' n )) 

»P, 


— <-—£' + —£ 


dt 2e 

where if we take eo, S small enough 

1 


0 < a = -(! - e 0 c' n - 6 ) < 1 f3 = 8\{ip + e 0 c' n ) . 


We combine with (15.101) . and thus consider 




whose solution is given by 
£ = 


So + fco 


1 — a 


£(zq) = £o 

P 


e — 


1 — a 


This curve provides a bound as long as d£/de < 0, i.e., for 


^max — 


i(l — a) (£{ 


1 — Ot \ 1 1 - 


1 — a 


< e < e 0 ;, 


which gives a maximal value of 

£max = £ (Ci 

(ignoring the dependence on v , A). 


G 


4— 2 a 
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This is somewhat unexpected. With respect to the scaling associated with the 3D 
NSE, the L 2 -condition is scaling-invariant, while the ^-Holder coherence condition 
exhibits sub-critical scaling. Hence, from this point of view, the scaling-invariant 
condition is a weaker condition. In contrast, when considering dynamics of the 
solutions, at least within the realm of bounding curves in the energy,enstrophy 
plane, the ^-Holder coherence condition appears to be much weaker as it allows for 
qualitatively higher maximum values of the enstrophy: exponential vs. algebraic (in 
Grashof). Incidentally, this is philosophically consistent with another very recent 
result illustrating discrepancies between the scaling and dynamical properties of 
solutions to the 3D NSE [5j. 

To compare the L? scaling-invariant case to the r-Holder coherence assumption 
in the case \ < r < 1, in terms of G, we set 


2r- 1 

and find that 
Since 


4r T 2 f f . eo Co T <5 -t- 5 

= 4 — 2a = 3 + eoCn + 1 5, i.e., r = 


2(eo c'q + S) + 2 
0 < eo Cq + 5 < 1 => 3/2<r<5/2. 

A r -L 2 

-> 4 — 2a for 1/2 < r < 3/2, 

2 r — 1 


we conclude that even the stronger ^ < r < 1 Holder coherence assumption-in 
principle-allows for larger values of the maximal enstrophy than the L^-integrability 
condition, and is in this sense a dynamically weaker condition. 
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